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Abstract 

The maximum intercluster distance and the 
maximum quantization error that are minimized 
by the MinMax algorithm are shown to be 
inappropriate error measures for color image 
quantization. A fast and effective (improves 
image quality) method for generalizing activity 
weighting to any histogram-based color 
quantization algorithm is presented. A new 
non-hierarchical color quantization technique 
called weighted MinMax that is a hybrid 
between the MinMax and Linde-Buzo-Gray 
(LBG) algorithms is also described. The 
weighted MinMax algorithm incorporates 
activity weighting and seeks to minimize 
WRMSE, whereby obtaining high quality 
quantized images with significantly less visual 
distortion than the MinMax algorithm. 

Key words and phrases: color image 
quantization, minimizing maximum intercluster 
distance, nonhierarchical clustering, image 
compression. 

1. INTRODUCTION 

In color quantization a truecolor image is 
irreversibly transformed into a color-mapped 
image consisting of K carefully selected 
representative colors. When Heckbert proposed 
the color quantization problem in his seminal 
SIGGRAPH paper [1], most graphic 
workstations had CRT monitors with an 8-bit 

O 

frame buffer, capable of displaying only K - 2 
= 256 colors at a time. With the advent of the 
Internet and multimedia eras, 24-bit frame 
buffers have become more commonplace; 


however, the same has made the need for 
efficient ways to represent images more 
important. Color quantization is a useful 
method for obtaining high-quality compressed 
images for storage and subsequent distribution 
via the World Wide Web. For a survey of some 
of the more popular color quantization 
techniques see [2, 3, 4, 5]. 

The remainder of the paper is organized as 
follows. Section 2 briefly describes the MinMax 
algorithm. A fast and effective method for 
generalizing activity weighting to any 
histogram-based color quantization algorithm is 
presented in Section 3. Section 4 presents a new 
non-hierarchical technique called weighted 
MinMax that incorporates activity weighing and 
is a hybrid between the MinMax and Linde- 
Buzo-Gray (LBG) algorithms. Section 5 shows 
that the quality of quantized images produced 
using the weighted MinMax algorithm is 
superior to the best version of Xiang’s MinMax 
algorithm. 

2. THE MINMAX ALGORITHM 1 

In 1985, Gonzalez [6] proposed an 
approximation algorithm that seeks to minimize 
the maximum intercluster distance. One of the 
interesting aspects of the MinMax algorithm is 
that Gonzalez [6] proved it constructs a AT-split 
whose maximum intercluser distance is less than 
or equal to two times optimal solution value. 

The application of Gonzalez’s general purpose 
MinMax clustering algorithm to color image 


* This algorithm is commonly referred to as either the MaxMin or 
MinMax algorithm; this paper shall use the latter. 
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Figure 1 : PSEUDOCODE FOR XIANG’S BEST VERSION OF THE MINMAX ALGORITHM. 

The shaded lines indicate changes from the original MinMax. 


quantization has been studied and enhanced by 
several researchers [7, 8, 9 and 10]. 

Pseudocode for the best version of Xiang’s 
MinMax algorithm [10] is provided in Figure 1. 
The basic steps of the original MinMax 
algorithm are as follows. Initially all of the 
colors in the image are in the same cluster, .v 0 , 
and its representative color, r 0 , is arbitrarily 
chosen. A new cluster is created in each of the 
K - 1 iterations of the primary for loop. A color 
that is furthest away from its representative color 
becomes the representative color of the new 
cluster. Every other color that is closer to the 
new cluster’s representative color than to its 
current representative color is moved to the new 
cluster. Xiang’s modifications to the original 
MinMax algorithm are summarized below: 

(i) the centroid of each of the K clusters 
produced in the for loop in lines 3-11 of 
Figure 1 become the representative colors for 
each of the K clusters (lines 12-13), and 

(ii) the distance function is scaled based on their 
relative contribution to luminance to partially 


compensate for the nonuniform nature of the 
RGB color space: 

d s (x,y) = 7(0 5 M - ,.R)J + (jt.G - y.Gj + (0 25 (i.fl - j.b)J . 

Neither of these modifications invalidates 
Gonzalez’s theorem that the AT-split produced by 
the MinMax algorithm has a maximum 
intercluser distance that is less than or equal to 
two times optimal solution value. 

3. SPATIAL ACTIVITY WEIGHTING 

The MinMax algorithm produces quantized 
images with considerable contouring because it 
treats all colors equally, whereas higher quality 
quantization techniques use a histogram to 
weight each color proportional to its frequency. 
In histogram-based, non-activity weighted 
algorithms, the importance of a color is based 
solely upon its frequency. That is, colors that 
occur more often in the image are given greater 
precedence. Various methods for weighting the 
importance of a pixel based on the spatial 



activity of the area of the image where the pixel 
occurs have been incorporated into color 
quantization algorithms [5, 11, 12, 13, 14, 15, 

16, 17]. This section describes a fast and 
effective method for generalizing activity 
weighting to any histogram-based color 
quantization algorithm [5]. 

The importance of a color is based upon its 
frequency and the activity level of the locations 
where the pixels occur by summing the activity 
weightings of pixels with the same color. The 
activity level of a pixel is calculated in a manner 
similar to Balasubramanian et al. [13]. The 
activity level a(x,y ) of a pixel at row x, column y 
of a truecolor image / is measured by the 
magnitude of the luminance gradient of the pixel 
given by: 

a(x,y) = I Y(I(x,y)) - Y(I(x-l,y)) I + 

I Y(I(x,y)) - Y(I(x,y+\)) I, (1) 

where Y is the luminance of the pixel, and 
assuming jc > 0 and y < H- 1 . When x is 0, the 
first term of Equation 1 is I Y(/(x, v)- Y(/(x+ 1 ,y ) I 
and when y isH-\, the second term of Equation 
1 is I Y(I(x,y)- Y(/(x,v-l ) I. The activity 
weighting w(x,y) of a pixel at row x, column y is 
calculated from its activity level using: 


H w : RGB R (3) 

where H w (c ) is the total activity weighting of the 
pixels in / with color c. By simply substituting 
the weighted histogram for the histogram, any 
histogram-based algorithm becomes activity 
weighted. Thus, a non-activity weighted, 
histogram-based algorithm may be viewed as an 
activity weighted, histogram-based algorithm in 
which all pixels are equally weighted as one. 
Activity weighting in this manner has been 
shown to be a fast and effective way to enhance 
the quality of images produced by histogram- 
based color image quantization techniques [5]. 

4. THE WEIGHTED MINMAX 
ALGORITHM 

The weighted MinMax algorithm uses the 
weighted histogram described in the previous 
section to incorporate activity weighting into the 
MinMax algorithm, whereby favoring colors 
with higher total activity weightings. 

Pseudocode for the proposed weighted MinMax 
algorithm is provided in Figure 2. The author’s 
modifications to Xiang’s best version of the 
MinMax algorithm given in Figure 1 are 
summarized below: 


1/4 

if a(x,y) = 0 

1/3 

if a(x> y) = 1 

!/a(jt, v) 

1 < a(x, v) < 12 

1/«U>) 125 

12 < a(x>y) < 16 

1/1 6 1 25 

if a(x , y) > 16 


(2) 


The author conducted a considerable amount 
of experimentation with various exponents and 
cut-off values to derive Equation 2. When the 
activity level is extremely low (0 or 1 ) the pixel 
is considered to be part of a flat region; 
therefore, it is considered less important than 
when its activity level is low (2). When the 
activity level is greater than 16, the pixel is 
considered to be part of an edge. Thus to 
prevent the blurring of edges, the activity 
weighting of a pixel is never less than 1 / 1 6 1 25 . 

The weighted histogram H w of a truecolor 
image / is defined as a total function: 


(i) the centroid is selected as the representative 
color for the initial cluster in line 2, in line 4 
the distance from each color, c, to its 
representative is weighted by H w (c), and 

(ii) the representative color of each cluster is set 
to its centroid at the end of each iteration of 
the primary for loop (lines 3-1 1) as in the 
LBG algorithm. 

Gonzalez’s theorem that the AT-split 
produced by the MinMax algorithm has a 
maximum intercluser distance that is less than or 
equal to two times optimal solution does not 
apply to the weighted MinMax algorithm for the 
following reasons: 

(i) the weighted distance used in line 4 of 
Figure 2 does not satisfy the triangle 
inequality, and 
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Figure 2: PSEUDOCODE FOR THE WEIGHTED MINMAX ALGORITHM. The shaded lines 
indicate changes from the Xiang’s version of the MinMax algorithm given in Figure 1. 


(ii) the cluster representatives may change at the 
end of each iteration of for loop in lines 3-1 1 
of Figure 2. 

5. EXPERIMENTAL RESULTS 

This section compares the effectiveness of the 
weighted MinMax algorithm to the MinMax 
algorithm in terms of both running time 2 and 
images quality. Twenty-five color images were 
used as an “in practice average-case” test set. 
The test set contains both natural and computer 
generated images, and contains a wide variety in 
both spatial and color resolution. In addition, 
the images in the test set exhibit many of the 
characteristics commonly found in color digital 
imagery: 

• low frequency (smoothly varying) regions 
where false contours are more noticeable, 


• isolated colors which force color shifting, 
and 

• high frequency regions where sharpness 
may be deteriorated. 

Image quality is difficult to quantify; 
however, the following distortion metrics are 
used: RMSE, WRMSE, MaxDiam, and 
MaxError. WRMSE weights the distance from 
each pixel in the quantized image to its 
corresponding pixel in the original image by the 
pixel’s activity weighting. The diameter of a 
cluster produced by either the MinMax or 
weighted MinMax algorithm is defined to be the 
maximum distance between any two colors in 
the cluster. MaxDiam is defined to be the 
maximum diameter of all K clusters. MaxError 
is the maximum distance between any pixel in 
the quantized image and its corresponding pixel 
in the original image. All of these error metrics 
use the Euclidean distance function, not the 
scaled distance function, d s . 


2 The algorithms were executed on a SUN Ultra I running Solaris 2.6 
with 256 MB of RAM when the machine was otherwise idle. 
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Table 1: EMPIRICAL DATA FOR MINMAX AND WEIGHTED MINMAX. K= 256, 24 bits of 
precision. The weighted MinMax data is shaded. 



(a) Time= 18.78s, RMSE=7.59, WRMSE=7.44, 
MaxError=30.56, MaxDiam^ 3.20 


(b) Tiroe= 106.65s, RMSE=6.85, WRMSE=6.57, 
MaxError=25S7, MaxDiam= 56.04 



(c) Time=29.43s, RMSE=5.38, WRMSE=3.58, 
MaxError=Sl .17 , MaxDiam=\21 .(>\ 


(d) Time=169.81s, RMSE=6.20, WRMSE=5.48, 
MaxError=4&.S5, MaxDiam=\ 12.54 


Figure 3: MINMAX (A-B) AND WEIGHTED MINMAX (C-D) QUANTIZATION OF SOUDS 
AND LENA. K= 256, 24 bits of precision 


Table 1 provides empirical data for both the 
MinMax and weighted MinMax algorithms. The 
average running time of the weighted MinMax 
algorithm is about 1 .6 times that of the MinMax 
algorithm. However, Figure 3 illustrates that 


the weighted MinMax algorithm produces 
higher quality quantized images than the 
MinMax algorithm. This is inconsistent with 
the observations that the MaxError and 
MaxDiam of the images quantized in Figure 3 


























using the weighted MinMax algorithm is larger 
than those quantized with the MinMax 
algorithm. Thus, neither the MaxError nor the 
MaxDiam metric is consistent with visual 
distortion. On the other hand, both the RMSE 
and the WRMSE of the images quantized in 
Figure 3 using the weighted MinMax algorithm 
is smaller than those quantized with the 
MinMax algorithm. Hence RMSE and 
WRMSE are more appropriate error measures 
than MaxError and MaxDiam are for color 
image quantization. 

6. CONCLUSION 

This paper conducted a comprehensive 
comparison between the MinMax algorithm and 
the proposed weighted MinMax algorithm. The 
comparison showed that MaxDiam and 
MaxError are not suitable distortion metrics for 
color image quantization. While the running 
time of the weighted MinMax algorithm is 
slower than the MinMax algorithm, the quality 
of quantized images produced was shown to be 
superior. A generalized method for applying 
activity weighting to any histogram-based color 
image quantization algorithm was also 
presented. The generalized activity weighting 
method was shown to be a fast and effective 
way to enhance the quality quantized images. 
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